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We examine the full nonlinear dynamics of closed FRW universes in the framework of D-branes 
formalism. Friedmann equations contain additional terms arising from the bulk-brane interaction 
that provide a concrete model for nonsingular bounces in the early phase of the universe. We 
construct nonsingular cosmological scenarios sourced with perfect fluids and a massive inflaton 
field which are past eternal, oscillory and may emerge into an inflationary phase due to nonlinear 
resonance mechanisms. Oscillatory behaviour becomes metastable when the system is driven into a 
resonance window of the parameter space of the models, with consequent break-up of KAM tori that 
trap the inflaton, leading the universe to the inflationary regime. A construction of the resonance 
chart of the models is made. Resonance windows are labeled by an integer n > 2, where n is related 
to the ratio of the frequencies in the scale factor/scalar field degrees of freedom. They are typically 
small compared to volume of the whole parameter space, and we examine the constraints imposed by 
nonlinear resonance in the physical domain of initial configurations so that inflation may be realized. 
We discuss the complex dynamics arising in this pre-inflationary stage, the structural stability of 
the resonance pattern and some of its possible imprints in the physics of inflation. We also approach 
the issue of initial configurations that are connected to a chaotic exit to inflation. Pure scalar 
field bouncing cosmologies are constructed. Contrary to models with perfect fluid components, the 
structure of the bouncing dynamics is highly sensitive to the initial amplitude and to the mass of the 
inflaton; dynamical potential barriers allowing for bounces appear as a new feature of the dynamics. 
We argue that if our actual Universe is a brane inflated by a parametric resonance mechanism 
triggered by the inflaton, some observable cosmological parameters should then have a signature of 
the particular resonance from which the brane inflated. 



PACS numbers: 

I. INTRODUCTION 

The issue of the initial conditions of our present Uni- 
verse is connected to the problem of the initial singular- 
ity and to the possible solutions adopted to circumvent 
this problem, which lie in the realm of a quantum the- 
ory of gravitation. In fact we may consider that the ini- 
tial conditions of our present expanding Universe were 
fixed when the early Universe emerged from a Planckian 
regime and started its classical evolution. However by 
evolving back the initial conditions using Einstein clas- 
sical equations the Universe is driven towards a singular 
point where the classical regime is no longer valid. This 
is an indication that classical General Relativity is not a 
complete theory and in this domain quantum processes 
must be taken into account. Therefore initial conditions 
from which our classical universe evolved should crucially 
depend on the version of quantum gravity theory adopted 
to describe the dynamics in the neighborhood of the clas- 
sical singularity. This implies that recent observational 
results in Cosmology could in principle guide us in nar- 
rowing the possibilities of choices. Inflationary cosmol- 
ogy, for instance, although a highly appealing theoretical 



paradigm, relies on assumptions about how the universe 
emerged from the cosmic singularity. In this vein models 
for a preinflationary phase, including quantum correc- 
tions and consistent with the inflationary paradigm, are 
important to be examined. 

Among several propositions to describe the dynamics 
in this preinflationary semiclassical domain are, for in- 
stance, quantum loop cosmology[l| and the string based 
formalism of D-branes both of them leading to correc- 
tions in Einstein's equations and encompassing General 
Relativity as a classical (low energy) limit. In the case of 
spatially homogeneous and isotropic cosmologies, the ba- 
sic resulting distinction between the two approaches lies 
in the corrections introduced in Friedmann's Hamiltonian 
constraint: quantum loop cosmology leads to corrections 
in the kinetic energy term of matter fields while bulk- 
brane corrections lead to extra potential energy terms. 
In both cases we may have bounces in the scale factor 
corresponding to the avoidance of a singularity in the 
models. 

In the present paper we adhere to the string based 
formalism of Z?-branes. In this scenario extra dimen- 
sions are introduced, the bulk space, and all the matter 
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in the universe would be trapped on a brane with three 
spatial dimensions; only gravitons would be allowed to 
leave the surface and move in the full bulkQ. At low 
energies General Relativity is recovered but at high en- 
ergies significant changes are introduced in the gravita- 
tional dynamics. Our main interest here is connected to 
the high energy /quantum corrections that are dominant 
in the neighborhood of the singularity, resulting in a re- 
pulsive force that avoids the singularity and leads the 
universe to undergo nonsingular bounces. An elegant 
geometrical derivation of braneworld dynamics embed- 
ded in 5-dim spacetimes may be found in Refs. JU-JU 
where both high-energy local corrections as well as non- 
local bulk corrections on a Friedmann-Robertson- Walker 
(FRW) brane are analyzed. Bouncing braneworld mod- 
els were constructed by Shtanov and Sahni[6] based upon 
a Randall-Sundrum type action with one extra timelike 
dimension. 

We here examine the full nonlinear dynamics of spa- 
tially closed FRW preinflationary braneworld models 
with a massive scalar field (the inflaton) and several non- 
interacting perfect fluids. The matter fields evolve on the 
brane, where high- energy /quantum gravity corrections 
due to the bulk are included and implement nonsingu- 
lar bounces. We have previously approached analogous 
models with a radiation fluid plus a scalar field in the 
form of small perturbations, where the brane corrections 
were due to the radiation fluid only or to a phantom-type 
fluid0. 

With the full nonlinear dynamics new possibilities for 
cosmological scenarios arise, as nonsingular oscillatory 
bouncing cosmologies sourced with a pure scalar field, or 
with a scalar field plus several perfect fluid components 
that allow to model the effect of dark matter together 
with barionic matter in the gravitational dynamics. Such 
nonsingular oscillatory solutions have the theoretical ad- 
vantage of avoiding the problem of initial conditions at 
past infinity occurring with one-single bounce solutions 
and, further, are favored by entropy considerations. We 
make a detailed examination of nonlinear parametric res- 
onance mechanisms that are present in the full dynamics 
and turn these bounded oscillatory solutions metastable 
allowing the model to emerge naturally into an inflation- 
ary phase. We consider the restriction such mechanisms 
impose on the initial configurations so that the models 
may realize inflation. 

We organize the paper as follows. In Section II we 
give a brief introduction to the framework of the brane- 
world formulation, making explicit the assumptions used 
in obtaining the dynamics, and derive the full dynami- 
cal equations of the models. In Section III we describe 
some basic structures (as critical points, invariant planes 
and attractors at infinity) that constitute the skeleton of 
phase space and allow to organize the dynamics in phase 
space. In section IV we restrict the matter content of the 
models to a massive scalar field plus dust and radiation, 
that constitute a minimal set of ingredients appropriate 
for a preinflationary model, and analyze the constraints 



on the parameters of the model so that the dynamics may 
allow for bounded oscillatory bouncing solutions. In Sec- 
tions V and VI we make a semi-analytical approach to 
nonlinear resonance phenomena in the models, that may 
turn the oscillatory bounded solutions into metastable 
ones, with an inflationary behavior, and we construct 
the resonance charts of the dynamics. In Section VII 
we treat the case of bouncing cosmologies sourced by a 
pure scalar field. Section VIII is devoted to the exam 
of the dynamics of initial condition sets connected to a 
chaotic exit to inflation. Section IX contains the conclu- 
sions and final discussions. Throughout the paper we use 
units such that h = c = 1. 



II. THE MODEL AND ITS DYNAMICS 

Our task here is to derive the full nonlinear dynamical 
equations of the models, (actually a 4-dim autonomous 
dynamical system with one first integral) and to ana- 
lyze structure of the associated phase space. In the 
framework of D-brane formalism, we consider a closed 
Friedmann-Robertson- Walker (FRW) metric on the 4- 
dim braneworld embedded in a 5-dim conformally flat 
bulk. The matter content of the models is constituted 
of a scalar field <j> plus several noninteracting perfect flu- 
ids, each with equation of state pi = otiPi. These matter 
fields are constrained to propagate on the brane only. 

We start by giving a brief introduction to brane-world 
theory, making explicit the specific assumptions used in 
obtaining the dynamics of the model. We rely on refer- 
ences 6] -[5], and our notation basically follows Q. Let 
us start with a 4-dim Lorentzian brane £ with metric g a b, 
embedded in a 5-dim conformally flat bulk A4 with met- 
ric gAB- Capital Latin indices range from to 4, small 
Latin indices range from to 3. We regard £ as a com- 
mon boundary of two pieces Aii and A^2 of A4 and the 
metric g a b induced on the brane by the metric of the two 
pieces should coincide although the extrinsic curvatures 
of a in M.\ and M.2 are allowed to be different. The 
action for the theory has the general form 

S = -L[ / ( < 5 )i?-2A 5 ) +2e / Ki 
+ f ( - 2A 5 ) - 2e / K 2 

JM.2 V ' •/£ J 

+ \ L {i)R ~ 2ct ) + L Li[9abi p ' ^ (1) 

In the above is the Ricci scalar of the Lorentzian 
5-dim metric gAB on M, A5 is the 5-dim cosmological 
constant and ^R is the scalar curvature of the induced 
metric g a b on S. The parameter a is denoted the brane 
tension, and K5 and K4 are Einstein constants in five- 
and four-dimensions, respectively. The unit vector n A 
normal to the boundary S has norm e. If e = — 1 the 
signature of the bulk space is (— , — , +, +, +), so that the 
extra dimension is timelike. The quantity K — K a b g ab 
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is the trace of the symmetric tensor of extrinsic curva- 
ture K ab — Y a c Y ih D V c n D , where Y A (x a ) are the 
embedding functions of £ in M. [6| • Also 



(2) 



is the Lagrangean density of the four dimensional mas- 
sive inflaton field <fr plus the perfect fluids (with equation 
of state pi = otiPi), whose dynamics is restricted to the 
brane S. They interact only with the induced metric g a b- 
We further assume that the inflaton field is nonminimally 
coupled with g ab , with coupling parameter £. All inte- 
grations over the bulk and the brane are taken with the 
natural volume elements y — e^g d 5 x and y— ( 4 )g ci 4 x 
respectively. 

Variations that leave the induced metric on X intact 
result in the equations 



(5) 



Gab + A 5 g A B = 0, 



(3) 



while considering arbitrary variations of gAB and taking 
into account ([3]) we obtain 

^G ab + e ^(sg ) -S^)= K ,(r ab -a 9ab ), (4) 

where S ab = Aafc — Kg ab . T ab is the energy-momentum 
tensor of the matter fields on the brane, resulting from 
@. In the limit k 4 — > oo equation (0J reduces to the 
Israel-Darmois junction condition^ 

( S ab - S ab) = 6 K 5 (r ab - CT 5ah ) . (5) 

We impose the ^-symmetryQ and use the junction con- 
ditions (|5|) to determine the extrinsic curvature on the 
brane, 
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(Tab - o T 5ab) 



(6) 



Now using Gauss equation 

{i) Rabcd = ^R M NRsY^Y^Y^Y s d + e(K ac K bd - 

KadKbcj together with equations © and ([5]) we arrive 
at the induced field equations on the brane 



to be considered. Accordingly Codazzi's equations imply 
that 



D a ri = 0, 



(9) 



where D a is the covariant derivative with respect to 
the induced metric g a b- Together with the contracted 
Bianchi's identities in ([7]) it results D a jr b = 0, which cor- 
responds to a sufficient condition for g a b to be embedded 
in a conformally flat bulk. 

Equations |7J and ([9]) are the dynamical equations of 
the gravitational field on the brane. In Eq. ([7} we iden- 
tify 



A, 



As 
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1 
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respectively as the effective cosmological constant, and 
the effective Newton's gravitational constant in the 
brane. Both depend basically on the brane tension a. 
Eq. ([7]) is similar to Einstein equations in 4 dimensions, 
except by the second term in the RHS that is a correc- 
tion resulting from the brane-bulk interaction quadratic 
in the extrinsic curvature, while Eq. ([9]) is the conserva- 
tion law for the matter on the brane. We recall that for 
the evaluation of the extrinsic curvature (|6]) we use the 
energy-momentum of the matter fields on the brane. In 
our model they are described by the Lagrangean density 
©• 

Let us consider the FRW metric on the brane given by 
the line element 



ds 2 = -dt 2 + a{tf 



dr 2 



1 - kr 2 



+ r 2 {d0 2 + sin 2 Odcj) 2 ) , (11) 



where a{t) is the scale factor. For the noninteracting 
perfect fluids considered in ([2]), each with equation of 
state pi — otiPi, we have then 



Pi = 



,3(l+cO ' 



(12) 



where Ei is a constant of motion. The components of the 
tensor ir ab (cf. flSJ)) are now expressed as 



7100 = Y2 { Pi + + A } 



(13) 



(4) G ob 



In the above 

1 c 

^ab — ^T ac Tjj 
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~ e cr T ab + e k 2 TT ab . (7) 
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J^ TT ab + g gabTcdT Cd - — g a b T 2 . (8) 



We remark the absence of the conformal tensor projec- 
tion in Eq. ([7]) since the 4-dim brane is embedded in a 
conformally flat bulk, which is the case of the FRW brane 
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Pi + p<j> + 



A } 2ffij + l{Y,p*+p* + A } x 



where 



Pi + 



2£ 



(m 2 + £i?)0 2 - <j) 2 + 



^ 2 2- + - + - 

a a A a A > 



9ij, 



4> z + m 2 4> 2 



P<l> 



(14) 



m 2 (j) 2 ), (15) 
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A = 3£ 



r/d 2 


k \ 
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h2-# 




a z / 




a J 



(16) 



a/a + (a/a) 2 + /c/a 2 



and the curvature scalar i? = 6 
The equations of motion (O and © reduce then to 

2 k A 4 8nGpf 



/ a \ z fc A 4 
W + ^~T 3 



A + ^ 2 + 



Pi + \$' 



-~m 
2 



2 i2 



---u-) 2 +4> 

a a z 



6^ 
a 



(17) 



■ + 3-<j) + m 2 (b + 6f 
a 



0, (18) 



respectively. Eq. ([T8|) is the Klein-Gordon (KG) equation 
for the inflaton field <fi while Eq. (fTTj) is the (0, 0) equation 
in ((Tj) . Eq. (fTT)) is in fact a first integral of the (z, j) 
equations provided the KG equation holds. 

The nonminimal coupling of the inflaton with gravi- 
tation considered here is partly motivated by quantum 
calculations in curved spacetimes (taking into account 
quantum backreaction, renormalization, etc.) and partly 
by enlarging the possibilities of constructing successful 
inflationary and preinflationary scenarios(cf. for instance 
[ic| ) . The case £ = is the usual minimal coupling of the 
scalar field with gravitation, and £ = 1/6 is the so-called 
conformal coupling 11]. 

It is not difficult to see from Eq. (|17l) that the choice 
e = — 1 (corresponding to the extra dimension with time- 
like character) will implement nonsingular bounces in the 
dynamics of the scale factor a, implying that the mod- 
els are nonsingular. In this case the brane tension a is 
required to be negative in order that the subsequent evo- 
lution of the universe be compatible with observations 

(cf. Eq. (EH)). 

In the remaining of the paper we will restrict our- 
selves to the case e = — 1 and positively curved FRW 
universes (fc > 0), as we are interested in preinflation- 
ary nonsigular dynamics with metastable oscillatory be- 
haviour. Partly for analytical and numerical simplicity 



we will also restrict our analysis to the case of conformal 
coupling £ = 1/6. The basic features of this case encom- 
pass that of the minimal coupling case £ = 0. Arbitrary 
£'s correspond to a higher dimensional parameter space, 
the scope of which is beyond the purpose of the present 
paper. The case of pure scalar field bouncing cosmolo- 
gies will demand a separate treatment, and will be dealt 
with in Section VII. Finally, for numerical computation 
purposes we fix n 2 = 6 and k = 1. 

Within the above restrictions, using the conformal 
time variable r (defined by dr = dt/a) and the rescaled 
scalar field ip = a<fi, Eqs. (p~7j) and ([T8j) simplify to 



a' 2 + a 2 
where 



a 4 4 kL, \ i 

-a =-T(p i , l p)-— A 



T( Pi , <p) = [a 4 Pi + \ (V 2 + (1 + ™ 2 « 2 K 



and 



tp" + (1 + m 2 a 2 ) (p = 0, 



(20) 



(21) 



where a prime denotes derivative with respect to the con- 
formal time t. Eqs. (fT9|) and (|21[) constitute the basic 
dynamical equations to be dealt with in this paper, de- 
fined in the phase space (a, ip, a', ip'). 

We note that for m = the dynamics is separable and, 
as a consequence, integrable. The Klein- Gordon equation 

(f2"Tj) has the first integral £° = {^p' 2 + (p 2 ^j /2 inthis case. 



III. THE SKELETON OF PHASE SPACE 



From the above equations (|T9l) and (|2T|) we derive the 
dynamical system 

<f' = Pip, 
p v ' = -(1 + to 2 a 2 ) tp, 
a = p a /6, 



Pa 



-6a + 4A 4 a 3 



KB 3 



< i ; - - I ) — tt~ m 2 aif 2 



(22) 



The first and third equations (|22p are mere redefinitions. 
(p a ,a) can be shown to be canonically conjugated (cf. 
description of the invariant plane dynamics, for instance). 



This is not the case of (p v , ip), since the first integral p^|) 
is not a Hamiltonian constraint in these variables due the 
presence of fourth-order powers in ip' . 
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Three basic structures organize the dynamics in the 
phase space of the above dynamical system, namely, (i) 
an invariant plane, (ii) critical points and (iii) separatri- 
ces, allowing us to give a global description of the motion 
of the models. 

The invariant plane is defined by 

tp = Q, Pip ip' = 0, (23) 

where the dynamics is integrable; orbits with initial con- 
ditions in this plane are totally contained in it, actually 
being similar to the dynamics in the sector (a,p a ) of the 
separable case. A first integral of motion is given by the 
Hamiltonian constraint 

^ + V(a) - |a|£ rad = (24) 

where the potential V(a) is defined 



The equations of motion are equivalent to Hamilton's 
equations generated from (|24p , corresponding to the third 
and fourth equations (|22[) restricted to the invariant 
plane ([2"5]). 

A careful examination shows that for perfect fluids 
with —1/3 < oi; L < 1 the last term in the potential ([23)1 
acts as an infinite potential barrier and is responsible 
for the avoidance of the singularity a — 0. These po- 
tential corrections are equivalent to fluids with negative 
energy densities. This is in accordance with the fact that 
indeed quantum effects can violate the classical energy 
conditions, and may avoid curvature singularities where 
classical general relativity breaks down[12j. Such quan- 
tum violations tend to occur on short scales and/or at 
high curvatures, which is the case in our present models. 

Critical points in the finite region of phase space are 
stationary solutions of Eqs. (|22[) . namely, the points 
(j?tp = 0,f = Q,Pa = 0,o = Ocrit)) at which the RHS of 
Eqs. ([22]) vanish. Obviously the critical points are con- 
tained in the invariant plane, with a = a cr it being the real 
positive roots of V'(a) = 0, where a prime here denotes 
derivative with respect to a. Depending on the values of 
the parameters ( A4, a, Ei, a,) we may have from one to 
several real positive a cr it- However it is not difficult to 
verify that not all of them satisfy the constraint equation 
(fT9|) and that at most two critical points are present, as- 
sociated with one minimum and one maximum of V(a). 
As a matter of fact, this is the case for a fixed A4 and 
sufficiently bounded values of Ei. These limiting condi- 
tions on the Ei have not in general a closed analytical 
form except for the cases of dust or radiation, and will 
be dealt with in next Section. For A4 = the dynami- 
cal system has one critical point only, corresponding to a 
minimum of the potential V(a); such configurations are 
not of interest for us since they corresponding to perpet- 
ually nonsingular oscillating universes, a scenario where 
inflation cannot be realized. 



To examine the nature of the critical points, and conse- 
quently the nature of the motion in their neighborhood, 
we linearize the dynamical system (|2"21 about each a cr it- 
The resulting 4x4 constant matrix of the linearization 
has the four eigenvalues 

A 1>2 = ±i^l + m^~ t , A 3 , 4 = ±^-V"{a crit )/6 (.26) 

related to the linearized motion about a crr t in the sec- 
tor ((f,p v ) and in the sector (a,p a ), respectively. This 
characterizes the minimum of V(a) as a center Pq, with 
two pairs of complex conjugated imaginary eigenvalues 
(F" (ac r it) > 0), and the maximum Pi as a saddle- 
center [l3j with one pair of real eigenvalues with opposite 
signs (V"(a cr it) < 0) and one pair of imaginary eigenval- 
ues. These results can be easily interpreted if we expand 
the integral of motion (fT9|) about the critical points as 

n = ^Pl + l V "( a ^ (« -«c r it) 2 

~ ^(^ + (l + m 2 a 2 rit V 2 ) 
+ E crit - \a\E rad + 0(3) = 0, (27) 

where 0(3) denotes higher order terms in the expansion 
and i? cr it = V(a a it) is the energy of the respective critical 
point. In a small neighborhood of the critical point these 
higher-order terms can be neglected and the motion is 
separable into the two sectors with approximate constant 
of motions 

E{a) = j^P 2 a + 7}V"(a clit ) (a - a C rit) 2 , 

E {v) = \{pI + (l + m 2 acr it V), (28) 

with E (a) - \<t\E {v) + Ecrit - W\E rad ~ and \E clit - 
\p\E ra d\ small. It is immediate to see that the sector asso- 
ciated with always correspond to rotational motion 
in the variables (ip,p v ) about the critical points, while 
the sector associated with E^ corresponds to either (i) 
rotational motion or (ii) hyperbolic motion in the vari- 
ables (a, p a ) about the respective critical point; namely, 
the minimum of the potential Pq corresponds to a center 
and the maximum Pi corresponds to a saddle-center, as 
mentioned before. Pq and Pi define, respectively, a sta- 
ble and an unstable static Einstein universe. The stable 
Einstein universe configuration has no classical analogue, 
arising from the dynamical balance between the perfect 
fluid energy content and the negative energy density con- 
nected to the high energy /quantum corrections in the 
potential. We are now ready to describe the topology 
of the motion about the critical points and its extension 
to a nonlinear neighborhood of these points. This will 
be done in next Section for the case of a preinflationary 
model containing cold dark matter (dust), dark energy 
(a positive cosmological constant) and radiation. 

Finally from the saddle-center critical point (when 
present) there emerges a structure of separatrices S con- 
tained in the invariant plane. One of them tends to a 



V{a) = ^-K^-\a\ £ -^r + i(E^r) 2 -(25) 
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DeSitter attractor at infinity, defining a escape of or- 
bits to the inflationary regime. In fact, a straightfor- 
ward analysis of the infinity of the phase space shows 
the presence of a pair of critical points in this region, 
one acting as an attractor (stable DeSitter configuration) 
and the other as a repeller (unstable DeSitter configura- 
tion) . The scale factor approaches the DeSitter attractor 
as a(r) ~ (Co — t)^ 1 for r — > Co, where r is the confor- 

mal time, or as a(t) ~ exp^ty^ A 4 /3^ ) for t — > oo, where 

t is the cosmic time. 

IV. PREINFLATIONARY MODEL WITH DARK 
MATTER, DARK ENERGY AND RADIATION 

To proceed with our analysis further we will consider 
braneworld models whose matter content is restricted to 



a dark matter component (dust) and radiation, together 
a massive scalar field (the inflaton) and a dark energy 
component described by the effective cosmological con- 
stant in the brane (cf. Eq. (jlO[))- The reasons for these 
restrictions are two-fold. First the model contains a min- 
imal set of ingredients that is appropriate for a preinfla- 
tionary model. The second is more technical and has to 
do with the number of independent perfect fluid compo- 
nents. A large number of fluid components results in a 
higher dimensional parameter space, what turns the nu- 
merical/analytical analysis of parametric resonance quite 
involved. Furthermore the resonance patterns that are 
important for the physics of inflation in this model are 
typical for the dynamics of the general model as we will 
discuss later. In this instance Eqs. (|19 p -(|2i p reduce then 
to 



ff + V(a) - |a|P rad - J£l (V 2 + (1 + m W 2 ) + ± + f** ) (f + (1 + m 2 a> 2 ) 

|i(/ + (l + W 2 aV) 2 = 0, (29) 



and 

if" + (1 + m 2 a 2 ) (f = 0, (30) 

where 

V{a) = 3a 2 - A 4 a 4 - \a\E dust a+U^- + .(31) 

2 V a A a J 

We remark that for the integrable case m = 0, Eq. (|29|) 
simplifies to 

EL + V(a)-\a\(E I3d + E° vl )=0, (32) 

where 

V(a) = 3a 2 - A 4 a 4 - \cr\E dust a 

+ \( E ^K ^)\ (33) 
2 V er a / 

with £® = ((p' 2 + <p 2 )/2 a constant of motion. If we 
compare Eq. (|32|) with Eq. (|24|) for the invariant plane, 
and Eq. (1551) with Eq. (|3"Tj) we may see that the inte- 
grable dynamics (m = 0) is analogous to the integrable 
dynamics in the invariant plane up to the substitution 
E ra d — > E'rad + £%■ In other words, the integrable scalar 
field behaves as a radiation fluid in respect to the dy- 
namics of the scale factor, a fact that will be used in 
Section VII, for the case of pure scalar field cosmology. 
We should remark that in the low energy limit the cos- 
mological constant on the brane A4 may be interpreted 



as the effective vacuum energy of the inflaton field and 
tp are the spatially homogeneous expectation values of 
the inflaton fluctuations about its vacuum state. When 
these fluctuations arc considered small or initially small, 
namely, taken near the invariant plane, we may neglect 
£° and the integrable dynamics m = is approximately 
that of the invariant plane. The scalar field fluctuations 
will then actually have the role of just triggering the res- 
onances in the perfect fluid cosmologies. 

In Fig. 1 we depict the phase space portrait of the 
integrable dynamics in the invariant plane ((p = = p v ) 
for the case of dust and radiation, with |cr|i?dust suffi- 
ciently bounded so that V(a) has a well (we adopted 
\a\ = 500 and E ra d — E dust = 10~ 3 ). The critical points 
Po (center) and Pi (saddle-center) correspond to stable 
and unstable Einstein universes. Typically the model al- 
lows for the presence of perpetually bouncing universes 
(periodic orbits) in region I of Fig. 1, associated with 
motion in the potential well V(a) about the stable Ein- 
stein universe configuration Po- These configurations are 
basically the ones for which \a\ P ra( j < V(a max ). Re- 
gion I, understood as a nonlinear neighborhood of Po, is 
bounded by the homoclinic separatrix S emerging from 
the saddle-center Pi[l4j. Orbits in Region II are cor- 
respond to one-bounce universes. A separatrix emerges 
from Pi towards the deSitter attractor at infinity, defin- 
ing a escape to inflation. 

Finally it is worth mentioning here some basic struc- 
tural differences between the integrable dynamics in the 
invariant plane and that of the integrable case (Eq. (|3"2")l 
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V(a) will have this well provided that 



rad 



< 



2187 1 



2048 A 4 3 ' 



(34) 



5 o 




'Escape to cfeSifi 



0.5 1 



FIG. 1: The phase portrait of invariant plane dynamics with 
the critical points Po (center) and Pi (saddle-center) corre- 
sponding to stable and unstable Einstein universes. The pe- 
riodic orbits of region I describe perpetually bouncing uni- 
verses. Orbits in Region II are solutions of one-bounce uni- 
verses. A separatrix 5 emerges from Pi towards the deSitter 
attractor defining a escape to inflation. 



and Eq. (|3"0|) for m = 0) whenever ip{0) and/or p v (0) are 
not zero. The phase space portrait in the plane (a,p a ) 
is similar to that of the invariant plane shown in Fig. 1, 
however P\ and Po are no longer critical points but pe- 
riodic orbits. The integrable dynamics is not restricted 
to this plane but, in the case of region I, evolve on tori 
that are the direct product of the closed curves in region 
I to periodic orbits in the sector (ip,p v ). In particular 
the product of the separatrix S that bounds the region 
I times a periodic orbit with initial conditions fixed by 
the constant £® is a cylinder. The orbits on this cylinder 
coalesce (for times r — > ±oo) to the periodic orbit at Pi 
with the same constant of motion. This cylinder is said 
to be homoclinic to the periodic orbit at Pi. 

As we mentioned already, for a fixed value of A4 
the potencial V(a) has two extrema (one minimum 
and one maximum) for suitably bounded values of 
(| a I, -Erad, -Edust), corresponding to a well in the poten- 
tial. For the case of pure radiation we can show that 



The two extrema of V{a) for this case are critical points 
of the dynamics since V(a) > for all a > imply- 
ing that even the minimum of the potential belongs to 
the physical phase space domain. For .E^ad violating the 
above restriction no extrema are present and the system 
has no critical point in the finite region of phase space. 
For the case of pure dust V(a) has no extrema if 



> 
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For radiation plus dust we do not have a closed form 
for the constraint conditions on the parameters. How- 
ever it is not difficult to check that, in general, the in- 
crease of |<r|-Edust, namely, the increase of the gravita- 
tional interaction strength of the dust component, has 
the effect of taking the minimum of V(a) out of the 
physical space leading eventually to a destruction of the 
well. In other words the increase of |<j|-Edust has the 
effect of reducing the phase space volume avaiable for 
bounded and/or initially bounded (metastable) configu- 
rations. This is illustrated in Fig. 2 where we plot the 
potential V(a) for several increasing values of |<x|Ed us t 
with fixed A4 and E ra d (as a matter of fact we fixed 
A 4 = 1.5, E dust = 0.001, E rad = 0.01 and varied |er|). We 
see that for |<r|-Edust ~ 2.3 the well in the potential is no 
longer present. 

The structure of motion about Po will be the main ob- 
ject of our interest in the following. This region is physi- 
cally more relevant than that of one bounce models since 
it avoids the theoretical problems of posing initial con- 
ditions at past infinity, this issue being possibly brought 
to the realm of a semi-classical quantization of the dy- 
namics in the reg ion I. Furthermore maximum entropy 
considerations [15j favor the stable Einstein universe as a 
suitable past configuration about which initial states of 
the prcinflationary universe oscillate. However we must 
provide a mechanism for such bounded and perpetually 
bouncing universe configurations to become metastable 
so that they may realize inflation and escape to the De- 
Sitter attractor at infinity. Such mechanism will be the 
nonlinear resonance provided by the interaction with the 
scalar field sector (ip,p v ) in the nonintegrable case, as we 
proceed to discuss. 

We should finally remark that the increase of |cr| -Edust, 
namely, the increase of Newton's constant on the brane 
and/or the increase of the dust content of the model, have 
the effect of taking the stable Einstein universe configu- 
ration out of the physical phase space (cf. Fig. 2); in the 
latter situation the physical phase space is not simply 
connected presenting a void with the structure of a solid 
torus. 
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is given by the third-kind elliptic integral]!*: 
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FIG. 2: Plots of the potential V(a) for increasing values of 
r\ = a | iSdust- For higher values of |o"|.Edust the minimum is 
out of the physical phase space and eventually disappears (for 
M-Edust ~ 2.3) together with the potential well. 



NONLINEAR RESONANCE OF KAM TORI 



We start by discussing the topology of motion about 
the stable Einstein universe Po that includes (i) the struc- 
ture of KAM tori that are present in the neighborhood of 
Po and arise from the tori of the integrable case to = 0; 
(ii) nonlinear resonance mechanisms inducing that either 
KAM tori are destroyed by the resonance allowing the 
orbit escape to the de Sitter attractor or the motion is 
stable (confined between two KAM tori) but otherwise 
stable. 

Let us consider the energy surfaces |cr|(P ra d + £") < 
Ed-it(Pi) = V(Pi) corresponding to bounded motion in 
the integrable case to = 0, or to initially bounded mo- 
tion in the nonintegrable cases. This phase space re- 
gion can indeed be characterized as a nonlinear neigh- 
borhood of the center Po, and is foliated by the two-tori 
S 1 x S 1 that are the topological product of periodic or- 
bits of the separable sectors (p>,p v ) and (a,p a ), with the 
two associated separately conserved quantities £ ° and 
£° a = (pg/12 + Via)) satisfying £° a - \a\(E rad + £°) = 0. 
The frequency v a of the periodic orbit in the sector (a,p a ) 



f>2 
03 



x dx 



(35) 



x* - 2a x x +(al + a*)) il b i=1 (x- fr) 



where (j3i,i = 1...6) and (ai izc^) are, respectively, the 
real and imaginary roots of V(a) — I^K-Erad + £") = 
( /3 3 < 02 < Pi are the positive real roots). The two tori 
of the integrable case are the topological product of the 
above class of periodic orbits parametrized by £° with the 
periodic orbits of the harmonic oscillator parametrized by 
£° (with frequency v v — l/2ir). 

For future reference we note that the periodic orbits of 
the sector (a,p a ) in the integrable case will be represented 
by the elliptic fixed point (tp — 0,^ = 0), namely the 
origin of the Poincare map with surface of section p a = 0, 
in case £ ° = 0. For £® ^ the corresponding integrable 
tori are represented by closed invariant curves about the 
origin of the map. For a small nonintegrable parame- 
ter to this picture is maintained with (tp = 0,p v = 0) 
as a center of a primary island of KAM tori; in fact the 
KAM theorem[l7| establishes the stability of tori with a 
sufficiently incomensurate frequency ratio, which in the 
present case means v a sufficiently irrational. Other in- 
tegrable tori are destroyed by the nonintegrable pertur- 
bation, and the region between two remaining invariant 
tori presents an intricate dynamics (unstable periodic or- 
bits, stable periodic orbits surrounded by islands, broken 
separatrices, and stochastic layers, this structure repeat- 
ing down to smaller scales[14j). The importance of KAM 
tori for Hamiltonian systems with two degrees of freedom 
comes from the fact that they prevent the diffusion of tra- 
jectories in the whole phase space, and thus preventing in 
our model the entrance to inflation. As m increases nu- 
merical experiments show that invariant KAM tori may 
be destroyed, with a consequent loss of stability of the 
system. This is the case of interest to us as orbits initially 
trapped about the center (tp — 0,p v — 0) can escape into 
the inflationary phase. An important mechanism for this 
break-up of invariant tori is nonlinear resonance that oc- 
curs in a restricted domain of parameters of the system as 
we proceed to examine using a semianalytical approach. 

If to is small and/or we start from an initial condition 
{potPipo) small (namely, near the invariant plane) we may 
approximate Eq. p0[) as 



p" + (l + TO 2 ag(T)) p = 0, 



(36) 



where ao(r) is a solution of a" + dV(a)/da « (cf. Eq. 
([22])). Eq. (|3"6"]) has the form of a Lame-type equation, 
and therefore parametric resonance occurs when the ratio 



R 



is a rational number, where lit! 



1 

2^ 



(0.9 to) 



{Pi 



PI) 



(37) 



(38) 
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Under this condition <p begins to grow exponentially in 
time and to act on the dynamics of the scale factor a 
which in turn will modify (|36|) . This feedback will re- 
structure the resonance, either (i) leading the dynam- 
ics into a more unstable behaviour, with amplification 
of the resonance mechanism and consequent breaking of 
the KAM tori and escaping of the orbits; or (ii) lead- 
ing the orbits to a general chaotic motion in a bounded 
region of phase space. This general nonlinear resonance 
mechanism can be given an approximate analytical treat- 
ment, through which we can fix the dominant reso- 
nances, namely, the rational numbers R corresponding to 
the dominant resonances. In this approximation, where 
higher-order terms in ((p,p v ) are neglected, the first in- 
tegral ([29]) is approximated by 



ME, 



rad • 



(39) 



Further, in the remaining nonintegrable term the vari- 
able a and p are substituted by the intcgrable solutions 
(Xq(t) and po(r). Since the Hamiltonian character of the 
first integral is recovered in the approximation, we are led 
to introduce action-angle variables {J a , ©a, Jip, 8b), the 



angle variables being defined as (9 Q 



such that both 6 a and <d v vary in the interval [0, 1] dur- 
ing a complete cycle of the original variables. We remark 
that in the numerical experiments considered we have 
£® ~ M-EVad in the initial stages of the dynamics. £® 
is however not conserved as the dynamics proceeds, the 
nonintegrable term being responsible for the exchange of 
energy with the sector (p,p v ). Taking into account that 
the function ao(r) is periodic with period T a — v^ 1 , the 
expansion in circular functions of the nonintegrable term 
of El) takes the formO 



Hi = — !|i m 2 J"i 0) j( 0) ( A n cos2n7re„) cos 4^, (40) 



where A n are numerical coefficients depending the pa- 
rameters of the model through the roots of V(a) — 
|c|-E ra d ~ (cf. Eq. ([35]) ). and a zero superscript de- 
notes action variables of the integrable case. Hamilton's 
equation for J a derived from (|39[) can be integrated for 
each given term n of the series, yielding to a first approx- 
imation 



3 a 



1^1 

2 



m 2 Ji 0) 4°' ^ 



A n 



cos(27mO a — 4ttQ v ) 



R-2/n 



+ 



cos(27rn0 a + 4-7r© ¥ 



R + 2/n 



(41) 



From (|41[) we have that the dominant terms are the 
ones for which R ~ 2/n. Therefore for a fixed n > 2 the 
expression 



R 



n > 2 



(42) 



determines a volume in the parameter space 
(\cr\, -Eradj Edust, m) in the neighborhood of which a 
n-resonance occurs. It represents a further step beyond 



the analysis of resonances (|37| in the linear regime of 
Lame- type equation (f3"u| . The set up of the resonance is 
signaled by the bifurcation of the periodic orbit at the 
origin {ip — 0,p v — 0) into an unstable periodic orbit 
plus one or two characteristic stable periodic orbits of 
the resonance (according respectively whether n is odd 
or even), a fact crucial to realize inflation, as we will 
discuss later. 

The resonance chart for the model can now be con- 
structed numerically using the exact dynamics. Expres- 
sion (|42[) - where approximations as well as the neglect of 
non-resonant terms were used - constitutes an accurate 
guide to localize and label the resonances in the parame- 
ter space (obviously for a fixed A4, the values of |er|£^dust 
and Eq must be compatible with initially bounded mo- 
tion). However, in the actual chart of resonance con- 
structed numerically using the exact dynamics, these do- 
mains will have a spread that is a correction of the ap- 
proximations due the full dynamics. 

The focus on the underlying bulk-brane structure of 
the gravitational dynamics will lead us to examine ini- 
tially the plane (a, m) of the parameter space, while fix- 
ing the total mass-energy of dust -Edust and of radiation 
EVad in several distinct ratios. For numerical purposes in 
all numerical experiments we will fix A4 = 1.5. 

Fig. 3 displays the resonant chart in the plane (a, m) of 
the parameter space for fixed -Edust = -Erad = 10~ 3 . The 
chart correspond to initial conditions taken near the in- 
variant plane, namely near the periodic orbit at the origin 
(ip = 0,p v = 0), with p a — p v — 0, tp — 1CP 4 . The con- 
tinuous lines are solutions of (l4"2l while the gray regions 
spreading about the lines are sections of the resonance 
windows by the planes E dust = 0.001, E Tad = 0.001. The 
remaining (white) regions are domains of stable motion, 
with the dynamics bounded by KAM tori. We restricted 
ourselves up to the resonance n = 5 in order to not over- 
crowd the Figure. On driving the system towards a res- 
onance zone (by an appropriate change of m and/or |er|) 
we can turn a stable configuration into a metastable with 
consequent escape of the orbits to the de Sitter attrac- 
tor at infinity in a finite time. In the realm of our pre- 
inflationary models, stability versus non-linear resonance 
instability will be considered connected to initial condi- 
tions near the invariant plane only, namely with (ip,p<p) 
small, corresponding to spatially homogeneous fluctua- 
tions of the inflaton field. Non-linear resonance will turn 
orbits generated from these initial conditions from stable 
to unstable (and vice versa) as a consequence of bifurca- 
tion of the critical point (ip — 0,p v — 0), at the origin 
of the Poincare map with surface of section p a = 0, from 
a center to a saddle (and vice versa), breaking up the 
KAM tori that trap the orbits about the origin. We re- 
call that the origin of the map is a periodic orbit of period 
l/u a in the (a,p a ) sector. This is illustrated by the two 
Poincare maps of Fig. 4 with surface of section p a = 0, 
for E = Ei = 10~ 3 and \a\ = 500 fixed, and m = 5.85 
in the domain of parametric resonance n — 3 (top map), 
and m — 7.1 in the domain of parametric stability be- 
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FIG. 3: Resonance chart in the plane (a, m) for -E ra d = 
Ed-ast = 10 -3 and ipo = 10~ 4 . The continuous lines are 
solutions of the approximate resonant condition (|42[) . while 
the gray regions about the lines are the parametric resonance 
windows for the exact dynamics, corresponding to the bifur- 
cation of the periodic orbit at the origin. The white regions 
correspond to KAM stable motion. 



tween the resonances n = 3 and n = 4 (bottom map). 
The origin in the top map is a saddle, with two associated 
centers, a consequence of the bifurcation of the periodic 
orbit at the origin due to the resonance n = 3. No KAM 
tori is present about the origin of the map so that orbits 
with initial conditions about the origin are not trapped 
and free to escape for large regions of phase space and 
eventually reach the deSitter attractor at infinity, realiz- 
ing inflation. The parametric domain of resonance thus 
favor inflation. The origin of the bottom map is a cen- 
ter, enclosed by KAM tori that trap the orbits with initial 
conditions in this neighborhood, forbidding escape to the 
deSitter infinity. Therefore the region of parametric sta- 
bility of the system will be unfavorable to the physics of 
inflation since the orbit (a configuration of the early uni- 
verse) will be trapped in a stable state between two KAM 
tori about the center at the origin. The structure of the 
stochastic sea is distinct in each case. If the system is 
in the region of parametric resonance, initial conditions 
near the invariant plane may undergo a long time diffu- 
sion through the stochastic sea to large regions of phase 
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FIG. 4: Poincare maps with surface of section p a = for 
(\a\ = 500, m — 5.85) in the domain of parametric resonance 
n = 3 (top), and (|<r| = 500, m = 7.1) in the domain of para- 
metric stability between the resonances n = 3 and n = 4 
(bottom). The origin in the top map is a saddle, connected 
to the bifurcation of the origin due to the resonance n = 3, 
favoring inflation. The bottom map has a center at the ori- 
gin, enclosed by KAM tori that trap the inflaton, preventing 
inflation. 
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FIG. 5: Resonance windows n — 3 for (1) dust plus radi- 
ation (_B rad = E dust = 1CT 3 ), (2) pure radiation (_E rad = 
10~ 3 ,S dust = 0) and (3) pure dust (-E rad = 0, E dust = HT 3 ). 
The range and relative size of the windows are distinct, the 
case of pure dust having a relative smaller volume though a 
more extended region of parameters. The presence of dust has 
the effect of squeezing the width of the resonance windows, 
reducing the domain of favorable configurations to realize in- 
flation as concerning the parametric resonance mechanism. 



space, and finally escape to the deSittcr attractor. On 
the other hand, when the system is in the region of para- 
metric stability initial conditions for orbits that diffuse 
are far from the origin, beyond the borders of the main 
island of the map, and diffusion with escape to deSitter 
infinity is extremely rapid. 

Numerical experiments with the exact dynamics show 
that the presence of dust has the effect of squeezing the 
volume of the resonance windows, being thus less favor- 
able to the occurence of inflation. In Figure 5 we show the 
resonance windows n = 3, corresponding to ipo = 10~ 4 , 
for fluid content with pure dust (-E ra d = 0, -Edust = 
0.001), pure radiation (E ra d = 0.001, .Edust = 0) and a 
mixture of both (-E ra d = -Edust = 0.001), the latter al- 
ready shown in Fig. 3. Their range in the parameter 
space and relative size are distinct, the avaiable width 
for a fixed a being much reduced when dust is present. 
For reference we give the width Am of n = 3 resonances 
windows for |cr| = 600 and |er|=400, and for several dis- 
tinct matter content: 



£ rad = 10" 3 £ rad = 10" 3 £ rad = 10" 3 £ rad = 0_ 
-Edust = -E dus t = 10~ 4 E dust = 10~ 3 E dust = 10 -3 

Am 1.1088 1.0446 0.51665 0.55880 

\a\ = 600 

Am 1.4418 1.3850 0.9307 1.1305 

| a | = 400 

Hence if we demand that dark matter is present as 
a dust fluid in this preinflationary phase then its total 
mass-energy content must not be exceedingly large in 
comparison with the radiation content in order that in- 
flation be properly realized. In addition we recall that, 
as discussed in Section IV, the presence of dust has the 
effect of reducing the avaiable phase space volume for 
bounded and/or initially bounded (metastable) configu- 
rations about the Einstein stable universe, as regulated 
by |cr|E dus t. 

Summarizing, the instability versus the stability of the 
origin (ip = 0,p v = 0) is crucial for the dynamics of in- 
flation, having a bearing on the dynamics of the spatially 
homogeneous expectation values <p(r) of the inflaton field 
related to the escape into inflation. In this instance the 
initial conditions for ip are assumed to be small, and are 
to be taken near the invariant plane (|23p. which cor- 
responds to a neighborhood of the critical point of the 
map at the origin (tp — = 0). Therefore the re- 

gion of parametric stability is unfavorable for producing 
inflation since the orbit (a configuration of the early uni- 
verse) will be trapped in a stable state enclosed by two 
invariant tori of a main KAM island of the map. On the 
other hand, on driving the system to a region of para- 
metric resonance, orbits with initial conditions near the 
invariant plane are turn into metastable configurations 
that either escape rapidly to de Sitter infinity or undergo 
a long time diffusion through stochastic regions of phase 
space before finally escaping. The resonance windows 
in the complete parameter space (tr, -E ra d, -Edust, ni) are 
constituted of n disjoint 4-dim volumes the sections of 
which - for instance with the surfaces -E ra d = 10 -3 and 
-Edust = 10~ 3 - result in the n gray regions of Fig. 3. 
The volumes of the windows are small as compared to 
the whole volume of the parameter space, and only ini- 
tial configurations inside them may realize inflation. We 
note that for a fixed m the larger the order n of the res- 
onance the stronger the gravitational interaction in the 
braneworld universe inflated from initial conditions con- 
nected with the resonance considered. 



VI. A PARTITION IN THE RESONACE 
WINDOWS CODED BY DISRUPTIVE 
RESONANCES 

From the point of view of the dynamics of inflation, 
the resonance windows present a further structure con- 
nected with disruptive resonances and/or long time diffu- 
sion before escape to inflation. In fact, as we proceed to 
discuss, a considerable domain of the resonance windows 
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- although corresponding to a bifurcation of the stable 
periodic orbit at the origin - does not lead to escape into 
inflation and must be properly discarded. To simplify 
our analysis we restrict ourselves to the resonance win- 
dow n = 3 (cf. Fig. 3) at fixed \a\ = 350. The associated 
range of m lies in the interval Am = [7.78431,8.84341]. 
A careful numerical examination shows that in this in- 
terval of the resonance zone we observe three dynam- 
ically distinct regions are present: (i) the domain on 
the left, ranging from the lower limit up to m < 8.52, 
corresponds to configurations for which the dynamics is 
highly unstable, the resonances being disruptive with a 
rapid escape to inflation (up to r < 15,000). (ii) A 
threshold region 8.52 < m < 8.56 that corresponds to 
configurations of orbits that undergo a long time diffu- 
sion (15, 000 < r < 100, 000) before escaping to inflation. 
(iii)The region on the right of the threshold, ranging from 
m > 8.56 up to the upper limit of Am. The motion of 
orbits connected to this latter parametric region - al- 
though corresponding to the case of a bifurcated saddle 
at the origin - is resonant and chaotic, but otherwise 
stable. Therefore region (iii) should be discarded since 
cosmological scenarios corresponding to these configura- 
tions do not properly realize inflation. 

The above dynamically distinct behaviors are illus- 
trated as follows. Fig. 6 shows the time signals for 
m = 7.99 taken in the region of disruptive resonance; the 
escape to inflation occurs at r ~ 140 and time signals are 
used since there is not enough recurrence for construct- 
ing a well defined Poincare map. Fig. 7 displays the 
Poincare map with surface of section p a = of a single 
orbit for m = 8.555, corresponding to the threshold re- 
gion (ii), undergoing a long time diffusion before exit to 
inflation, at r ~ 64, 000. Finally the time signals shown 
in Fig. 8 correspond to m = 8.6 in region (iii) near the 
right border of the resonance window n = 3. The motion 
is resonant and chaotic, but otherwise stable. All Figures 
were generated with ip' = = p a , <p = 10~ 4 . 

The above substructure appears to be a feature of the 
whole resonance window, as we have checked for other 
values of \a\ in the n = 3 window, as well as for other 
resonance windows. It is worth mentioning that in the 
case of pure dust the subdomain of disruptive resonances 
is a very thin sheet on the left border of the resonance 
window, while the larger portion of the window that al- 
lows for escape corresponds to long time diffusion. 

Finally we should remark that - as concerning nonlin- 
ear resonance phenomena - the general picture is that 
bouncing oscillating braneworld models have a small re- 
stricted domain in their parameter space where inflation 
can be realized. Typically variation of the parameters can 
shrink or stretch the resonance zones. However the un- 
derlying pattern of resonance windows and their internal 
substructure is maintained as we have checked numeri- 
cally. In this sense the pattern is said to be structurally 
stable. 
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FIG. 6: Time signals for m = 7.99 and \a\ — 350 correspond- 
ing to the region of disruptive resonances, near the left border 
of the n = 3 resonance window of Fig. 3. 
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FIG. 7: Poincare map with surface of section a' = of a 
single orbit, for m = 8.555 and \o\ — 350 in the threshold 
region of the n — 3 resonance window of Fig. 3. The orbit 
undergoes a long time diffusion before escape to inflation at 
r ~ 64, 000. The map exhibits the structure of the random 
motion of the orbit in the stochastic sea about primary and 
secondary KAM islands of the resonance. 



VII. PURE SCALAR FIELD BOUNCING 
COSMOLOGIES. METASTABLE DYNAMICAL 
CONFINEMENT 

We now proceed to examine the case when the model 
has no perfect fluid component. As we will see, even then 
the bulk-brane corrections allow for the presence of either 
perpetually oscillatory or metastable oscillatory models 
that emerge into an inflationary phase after a finite time. 
The dynamical equations in this case are given by Eq. 
m . namely the KG equation, 



if" + (l+m 2 a 2 ) <p = 0, 



S- o 




plus the modified Eq. JT 



where 



A 



4 4 

a 



AO) 



1 

6a 1 



A 2 M 



[ %)=0, (43) 



tp' + (1 + m 2 a 2 )tp 



(44) 



FIG. 8: Time signals for m = 8.6 and |a| = 350, correspond- 
ing to the region (iii) near the right border of the n = 3 
resonance window of Fig. 3. The motion is resonant and 
chaotic, but otherwise stable, with no exit to inflation. 
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For the integrable case m = 0, we have that Ao = 
W 2 (t) + <^ 2 (t)]/2 = [y/ 2 (0) + <^ 2 (0)]/2 is a constant of 
motion and Eq. (|4"3")l reduces then to 

a' 2 + V(a)= l ^-A a , V(a)=a*-^a i + ^L, (45) 

so that the dynamics of the scale factor a(r) is analogous 
to that of a radiation dominated bouncing braneworld 
universe with total energy content A (cf. Eqs. (TT9")) 
and (|2"T1) for _E rad ^ 0, other £^=0). Physical motion 
imposes the constraint V(a m - m ) < |<r|Ao, where T^(a m i n ) 
is the minimum of the potential (1451 , the inequality cor- 
responding to oscillatory motion. A distinct feature of 
the phase space is that the infinite barrier avoiding the 
singularity as well as the constant of motion |tr| Ao are 
both built up with the nonzero initial amplitudes of the 
inflaton field. A portrait of the phase space plane (a', a) - 
analogous to that in Fig. 2 - may be constructed by vary- 
ing |cr|A . In the domain V(a min ) < \a\A < V(a max ) 
the dynamics lies on tori that are the product of the peri- 
odic orbits in this phase plane times the circles ip' 2 +ip 2 = 
2Aq. The minimum configuration is a stable Einstein 
universe sourced by a massless inflaton. This stable Ein- 
stein universe is a limiting 1-dim torus, being actually 
the product of the point (a' = 0, a = a m ; n ) times the cir- 
cle defined by A = ^(a m i n )/|er|. Contrary to the cases 
with a perfect fluid component, the Einstein universe is 
not a critical point of the dynamics. No invariant plane 
is present. 

The compact phase space domain containing tori is 
analogously restricted by with the obvious substi- 
tution of Eg by Ao- We note that, for fixed a and A4, 
each of the integrable tori is generated by just one curve 
with initial conditions (ip' ,ipo). In this sense the energy 
surfaces |cr| Ao=const. contain just one curve. 

Analogous portrait can be constructed for a single orbit 
with initial conditions (^'g,^) by varying er, character- 
izing thus universes with distinct gravitational strengths 
Gn = I cr I / 87T. The corresponding tori yield qualitatively 
the same stable dynamics, about the stable Einstein uni- 
verse sourced by a massless conformal inflaton with a 
minimal gravitational strength. 

For the nonintegrable case, with m small (and 
\a\m 2 tp(0) 2 /6 sufficiently smaller than 1, as we discuss 
below), the pattern is that of stable motion on KAM 
invariant tori or between two invariant KAM tori re- 
sulting from the integrable case. The stability of mo- 
tion implies that all orbits are trapped between invari- 
ant tori and cannot escape to inflation, namely, infla- 
tion cannot be realized for these configurations. This is 
illustrated in Fig. 9 where we construct the Poincare 
map with surface of section a' = for fixed parameters 
(A4 = 1.5, \a\ = 500,m = 8.15). The innermost circle 
about the origin of the map is the section of the 1-dim 
torus (namely, the periodic orbit in the plane (<p',(fi) ) 
corresponding to the stable Einstein universe. In this 
example the stable Einstein universe - sourced by a mas- 
sive (m = 8.15) conformal inflaton - is described by 
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FIG. 9: Poincare map with surface of section a = for pure 
scalar filed cosmologies, with |oj = 500 and m = 8.15. The 
innermost circle about the origin is a 1-dim torus correspond- 
ing to a stable Einstein universe sourced by a scalar field. The 
domain inside this circle is unphysical. The last KAM torus 
shown has initial conditions (ip' — 0, ipo ~ 0.01105). Escape 
to inflation may occur in the outer border (for ipo > 0.01105) 
with oscillatory motion due to partial confinement of orbits. 



a' = 0, a ~ 0.00519690052) times the circle in the (<p', (p)- 
plane generated, for instance, from the initial conditions 
(ip'(0) = 0, if ~ 0.000696503226) with associated energy 
constant A ~ 0.2429935 x 10~ 6 . Below this scalar field 
energy no oscillatory motion is found, resulting in the 
void (no physical motion) about the origin of the map 
bounded by the above described innermost circle. The 
pattern of KAM tori sections, with eventual bifurcated 
secondary islands, extends up to initial conditions a' = 0, 
ip' = and ipo ~ 0.01105 where approximately the last 
confining KAM torus lies. In this region of trapped orbits 
we may note, for instance, 8 secondary islands (the far 
right centered about ipo = 0.01079) connected to a 8/5 
bifurcation. The outside region, denoted the stochastic 
sea in border of the main KAM island, corresponds to 
the domain of initial conditions where inflation can be 
realized in pure scalar field sourced braneworld cosmolo- 
gies. In this stochastic sea we may have either a rapid 
escape to inflation (for ipo — 0.0118 with r ~ 210, or 
for ipo — 0.01125 with r ~ 500, for instance) or a diffu- 
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sion with escape, for instance about the border of the two 
sets of 3 secondary islands connected to a 3/2 bifurcation 
(for ip = ±0.011305 with r ~ 1690). For larger initial 
conditions, namely (fo > 0.01234, no bifurcated islands 
are found in the stochastic sea but finite time oscillations 
are still found due to a mechanism of dynamical partial 
confinement, connected to values of the parameter 

<r = \a\m 2 (p(0f/6 (46) 

sufficiently close to 1, as we describe in the following. 

Let us recall that the potential V(a) in (|45|) presents 
a well (with possible oscillatory motion) for properly 
bounded values of A restricted by (|34p . Furthermore 
V(a) — |cr| Ao/3 = must have three real positive roots 
such that oscillatory motion is present. However the re- 
striction (f3~4"|) was derived for the integrable case and for 
the spatial curvature parameter rescaled to k = 1. Ac- 
tually the presence of the physical well is due to the bal- 
ance between the spatial curvature term and the infinite 
barrier term originated from the bulk-brane correction, 
more specifically, the balance between the values of the 
parameters k and Ao. It is not difficult to see that, in the 
integrable case, the decrease (increase) of k with Ao fixed 
as well as the increase (decrease) of Ao for k fixed, may 
destroy (or create) a well. Now for the nonintegrable case 
(to 7^ 0) a careful examination of the integral constraint 
(|4"3"|) shows that the term proportional to to 2 in the RHS 
will contribute to correct the curvature term to an effec- 
tive spatial curvature (fc e ff — 1 — <r) at the initial times. 
For the parameter configuration of Fig. 9 numerical eval- 
uations show that the above effects will be crucial for the 
dynamics of the scale factor a(r) when k c g < 0.157, that 
corresponds to <po •> 0.01234. In this instance, two possi- 
bilities arise: (i) for 0.01234 < ip Q < 0.1334, tp' = 0, the 
function 

p (a) ^ a 2_^ a 4 + _L A 2 M _M A( ^ )5 (47) 

(cf. (1431) ) has two extrema (one maximum and one 
minimum) but P(a) = has only one real positive 
root ao, meaning that the value of the relative maxi- 
mum of V(a) is smaller than zero. The orbit evolved 
from this initial condition would in principle be expected 
to escape without bounce. However, due to the in- 
crease of the effective time-dependent spatial curvature 
k c g ~ 1 — |(t|to 2 <^(t) 2 /6, the relative maximum of V(a) is 
raised above zero at a later time so that the orbit bounces 
back. This is illustrated in Fig. 10(a) where we plot 
■p(a) for the initial time with (ip' = 0,(fio = 0.0126) 
which has the only real root ao — 0.016955927198 (con- 
tinuous curve), and for the time of the first bounce of 
the orbit generated from these initial conditions (dashed 
curve) . The first bounce corresponds to the point (ip' ~ 
-0.0144932,(^0 ^ -0.00419929)) with k e g ~ 0.9. We 
note that the dashed curve represents a dynamical po- 
tential with its maximum above zero, allowing for the 
bounce. The adjustment of this effect with the period 
of the massive scalar field allows for a series of bounces 



of the orbit before it finally escapes, as shown in Fig. 
10(b). (ii) For ip > 0.1334 V(a) has no extremum but 
due to the effective time-dependent curvature two ex- 
trema are dynamically created, with the relative maxi- 
mum above zero, so that the orbit undergoes bounces in 
a relatively smaller number than in case (i) before es- 
caping. It is worth mentioning that this effect may be 
present for ipo < 0.01234, provided we properly increase 
the parameters to and/or |cr| to turn c; sufficiently small. 

Finally we should remark that the dynamics of pure 
scalar field brane cosmologies is not connected to a para- 
metric resonance pattern. The features of initial condi- 
tion domains that realize inflation in this case are similar 
to the ones for the domain of parametric stability (com- 
pare the Poincare maps in Fig. 4(bottom) and Fig. 9), 
as the configurations that realize inflation are the ones at 
the border of the primary KAM islands, corresponding 
to large values of ip. 



VIII. MOTION ABOUT THE 
SADDLE-CENTER. HOMOCLINIC CHAOS AND 
THE CHAOTIC EXIT TO INFLATION 

The saddle-center critical Pi induces, in the phase 
space of the models, the topology of stable and unstable 
homoclinic cylinders which emanate from unstable peri- 
odic orbits that exist in their neighborhood. This struc- 
ture was examined in some detail in [20l | and is briefly 
discussed here for completeness. The nonintegrability of 
the dynamics induces the breaking and crossing of the 
homoclinic cylinders leading to a chaotic exit to inflation 
as we proceed to show in this Section, where for simplic- 
ity we restrict ourselves to two perfect fluid components 
(dust and radiation only) as in Section IV. 

Our starting point is a fudamental property of saddle- 
centers given by Moser's theorem [21j. Moser's result 
states that it is always possible to find a set of canonical 
variables such that in a small neighborhood of a saddle- 
center the energy integral of motion (|19p is separable into 
rotational and hyperbolic motion pieces. The variables 
(a, p a , <p,p v ) of our system are already of the Moser type 
and, in a linear neighborhood of P±, we have the sep- 
arable rotational motion energy piece Et v ) and the hy- 
perbolic motion energy piece Er a ) as expressed in (j2"B")) . 
We recall that in this case V"(a cr i t ) < 0. In this ap- 
proximation, we note that the scale factor a(r) has pure 
hyperbolic motion and is completely decoupled from the 
pure rotational motion of the inflaton fluctuation tp. Let 
us consider the possible motions in this neighborhood. 
In the case E a = and p a = = (a — a cr it), the mo- 
tion corresponds to unstable periodic orbits T\ a \E rad m 
the plane (tp,p v ). Such orbits depend continuously on 
the parameter |<j|-E ra d. For E a = there is still the pos- 
sibility p a — ±\/6| V"(a cr it)| (a — a cr j t ), which defines the 
linear stable Vs and unstable Vjj 1-dim manifolds. The 
direct product of the periodic orbit T\ a \E lad with Vs and 
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FIG. 10: Plot of the function V(a) for the initial time (contin- 
uous curve) and the time of the first bounce (dashed curve). 
The occurrence of the bounce is due to the dynamical raising 
of the potential as a consequence of the increase of the effec- 
tive spatial curvature fc e ff(i~) as ip decreases. The first bounce 
occurs when k c a ~ 0.9. A sequence of this process produces 
the oscillatory behaviour shown below. 



Vu generates, in the linear neighborhood of the saddle- 
center, the structure of pairs of stable (7J CT |£; rild x Vs) and 
unstable (7| CT |E rad x Vu) 2-dim cylinders. We note that 
for -E cr it ~ l^i-Erad the manifolds Vs and Vjj are tangent 
at Pi to the separatrices S of the invariant plane con- 
taining the saddle-center; hence a pair of cylinders (one 
stable and one unstable) emanates from the neighbor- 
hood of Pi towards a ~ 0, while another pair emanates 
towards the two deSitter attractors at infinity (cf. Fig. 
1). Orbits on the cylinders coalesce into the periodic or- 
bit T\ a \E Tad asymptotically, the orbits being contained in 
the same energy surface |cr|B ra d as that of the periodic 
orbit; these orbits are denoted homoclinic to the periodic 
orbit 7[ CT |B rad or simply homoclinic orbits. Noticing that 



in general E( a ) ~ MB, 



\ V ) + B crit - \a\E rad ~ (cf. 423) 
and (|28|) and that E v is strictly positive we must have 
B cr it — | c | -Brad > 0; hence only energy surfaces satisfying 
this condition contain homoclinic cylinders. 

The nonlinear extension of the plane of rotational mo- 
tion (when 0(3) terms in ([27]) are taken into account) is 
a 2-dim manifold, the center manifold of unstable peri- 
odic orbits of the system parametrized with the energy 
I c I B ra d. The intersection of the center manifold with the 
energy surface |cr|B ra d is a periodic orbit 7T CT |E rad from 
which two pairs of cylinders emanate, as in the linear case 
previously described. It can be shown that in general the 
nonlinear extension of the center manifold folds, so that 
the unstable periodic orbits are no longer contained in 
the plane (ip,p v ). Now the extension of the cylinders 
away from the periodic orbits depends crucially on the 
integrability or nonintegrability of the system. 

In the integrable case (m = 0) the homoclinic cylin- 
ders are the topological product of the separatrices S 
times the periodic orbit with energy We are basi- 
cally interested in the pair of cylinders (one stable and 
one unstable) that emerge from the periodic orbit to- 
wards decreasing values of a, namely, in region I. It is not 
difficult to see that, in the integrable case, the unstable 
cylinder coalesces smoothly into the stable one. How- 
ever when the nonintegrability is switched on (namely, 
in 7^ 0), the smooth continuation of the unstable cylin- 
der into the stable one breaks, inducing the transversal 
crossing of them. The points of intersection define homo- 
clinic orbits which emerge from the periodic orbit along 
the unstable cylinder and return to it along the stable 
one, in an infinite time. Typically if the 2-dim cylinders 
intersect transversally once, they will intersect each other 
an infinite number of times, producing an infinite set of 
homoclinic orbits that are bi-asymptotic to the unstable 
periodic orbit of the center manifold. This infinite set 
of homoclinic orbits is denoted the intersection manifold. 
The dynamics near homoclinic orbits is very complex, as- 
sociated with the presence of the well-known horseshoe 
structures (cf. [lj|, [13] and references therein), with 
the homoclinic intersection manifold giving origin to the 
homoclinic tangle which is a signature of chaos in the 
model. 

This manifestation of homoclinic chaos is illustrated in 
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Fig. 11 where the fractality of the initial conditions basin 
boundaries, connected to the code recollapse/escape into 
inflation is shown. We select the initial condition set 
whose projection on (tp,p,p) is a square of characteris- 
tic length R ~ 10 -3 , constructed about the point (a = 
0.15, p a = 3.114327925239397, <p = 0,p v = 0) closed! to 
the separatrix S reaching Pi, containing N — 160,000 
initial conditions. These points obviously satisfy the 
constraint Eq. (j29l) . with Pdust = -E*ad = 10 ~ 3 and 
\a\ = 762, such that \<j\E rad is slightly below P crit = 
Knax- We note that these sets correspond to initial con- 
ditions for expanding universes that visit a neighborhood 
of the saddle-center Pi before either to recollapse or to 
escape into inflation. The points of the square are color- 
coded according to their asymptotic behaviour for a time 
r = 4, 000 - white corresponding to escape into infla- 
tion and black corresponding to recollapse to the neigh- 
borhood a ~ - resulting in the plots of Fig. 11 for 
m = 6 and m = 18. We notice that as we increase m, 
namely, as the system becomes more nonintegrable, the 
dominance of escape into inflation and fractality of the 
boundaries increase, characterizing a chaotic exit to in- 
flation. The dynamics involves no amplification of the 
inflaton field (namely ip and/or p v ), as opposed to the 
metastable nonlinear resonance regime before escape to 
inflation, described in Sections V-VI. 

The entanglement and infinite transversal crossings of 
the stable with the unstable homoclinic cylinders engen- 
der a further mechanism in the chaotic exit to inflation 
which we denote as draining of initial condition basins. 
Let us remark on the fact that the surface of the cylinders 
constitute a boundary for the general flow. If we start, 
for instance, with a set of initial conditions corresponding 
to initially expanding universes, two distinct flows will 
be associated with these initial conditions depending on 
whether they are contained inside or outside the stable 
cylinder. The flow corresponding to initial conditions in- 
side the stable cylinder will reach a neighborhood of the 
saddle-center Pi (with E a < 0) and will return towards 
the neighborhood of a ~ inside the unstable cylinder, 
while the flow of orbits associated with initial conditions 
that are outside the cylinder will reach the neighborhood 
of Pi and escape towards the deSitter attractor along the 
exterior of the unstable cylinder of the second pair. Now 
consider the first transversal intersection of the cylinders: 
a part of the orbits inside the unstable cylinder will enter 
the interior of the stable cylinder and the flow will pro- 
ceed inside the stable cylinder towards the neighborhood 
of Pi, from where it will reenter the unstable tube and 
proceeds towards the region a ~ and by a new inter- 
section a part of these orbits will again enter the stable 
tube and proceeds back towards the neighborhood of Pi 
and so on. The portion of orbits that in this subsequent 
intersection remained outside the stable cylinder will also 
proceed along it towards Pi and will escape to inflation. 
As the motion proceeds, the successive infinite intersec- 
tions of the cylinders drain the basin of initial conditions 
from black to white as r increases, as illustrated in Fig. 
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FIG. 11: Fractal basin boundaries, in the plane (p,p v ), of ini- 
tial condition sets for initially expanding universes that visit 
a small neighborhood of the saddle center before either recol- 
lapsing or escaping to inflation, in the case Eg = Ei = 10 -3 . 
Black points correspond to orbits that recollapse to the do- 
main a ~ and white points correspond to orbits that escape 
to inflation, for a time r = 4, 000, and m = 6 (left) and 
m = 18 (right). 



12, where we start from the initial condition set of Fig. 
11 for m = 6, color-coded at r = 4,000 (in the present 
Figures we used N = 40, 000 initial conditions only). The 
only orbits remaining in an infinite recurrence of the mo- 
tion are the homoclinic orbits, which constitute the ho- 
moclinic intersection manifold, and the infinite countable 
set of periodic orbits with arbitrarily long periods in the 
neighborhood of each homoclinic orbit. These orbits con- 
stitute a Cantor set that is a topological characterization 
of chaos in the model [l4|. 

Finally in Fig. 13 we exhibit two initial condition sets 
for m — 6, color-coded at r = 4, 000, in the cases of 
pure dust (left) and pure radiation (right), both with 
N = 160,000 points. For the radiation case the ini- 
tial set was constructed about the point (do = 0.15, 
p aa = 4.144304603687830, p = = p v0 ), sufficiently 
close to the separatrix, with \a\ = 1499 and P ra d = 10~ 3 . 
For dust the initial set was constructed about the point 
(a = 0.15, p a = 1.460426592928701, ip = = p v ), also 
sufficiently close[23j to the separatrix, with \a\ = 1630 
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FIG. 12: Illustration of the process of draining (in time) the 
initial condition basin of trapped orbits (black) in favor of 
escaping orbits (white), for m = 6. The initial condition set 
taken about the separatrix is the same as in Fig. 11. The 
Figures show the almost complete exit to inflation (up to a 
Cantor set of orbits) in a long time term. 



FIG. 13: Fractal basin boundaries in the plane {f,p v ) for 
pure dust (left) and pure radiation (right), for m = 6 and 
color-coded at r = 4, 000. The fractalization of the bound- 
aries appears less effective in both cases than in the case of 
radiation plus dust of Fig. ll(left). 



and -Edust = 10" 3 , both sets also with characteristic 
length R ~ 10~ 3 . Thus the mechanism of chaotic exit 
to inflation is enhanced in the case of dust plus radiation 
(see Fig. 11), in comparison with the case of pure dust 
or pure radiation components. 

The chaotic exit to inflation examined in the present 
Section is not a feature of bouncing inflationary cos- 
mologies only. In fact a similar pattern for fractal ini- 
tial conditions basin boundaries connected to the escape 
to inflation was observed and discussed by Cornish and 
Levin [24| in the context of singular closed FRW cosmolo- 
gies sourced with several conformally and/or minimally 
coupled scalar fields. An extensive analysis showed that 
the pattern is typical for the nonintegrable cases of the 
models, and a quantitative measure of the fractality was 
made by evaluating the box-counting dimension of the 
boundaries. The effect of draining of initial conditions 
from recollapse to escape as time increases, observed 
in our models, is not seen there due to the absence of 
bounces that avoid the singularity and allow for a long 
time recurrence. The Cantor set resulting from this effect 
as r — > oo corresponds to the strange repeller structure 



already appearing in their dynamics. 



IX. CONCLUSIONS AND FINAL DISCUSSIONS 

In the present paper we construct nonsingular cos- 
mological scenarios, in the realm of string inspired 
braneworld models, which are past-eternal, oscillating 
and may emerge into an inflationary phase due to non- 
linear resonance mechanisms. We consider a closed 
Friedmann-Robertson- Walker (FRW) metric on the 4- 
dim braneword embedded in a 5-dim conformally flat 
bulk. Local bulk effects on the 4-dim FRW braneworld 
introduce corrections in Friedmann's equations that al- 
low to implement nonsingular bounces in the scale factor 
of the models. This is the case of a timelike extra di- 
mension, when the corrections result in a repulsive force 
that avoids the singularity and provide a concrete model 
for bounces in the early phase of the universe. The mat- 
ter content of the model, confined to the FRW brane, 
consists of noninteracting perfect fluids with equation of 
state pi = ai Pi, with — 1/3 < on < 1, plus a massive con- 
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formally coupled scalar field. In the FRW brane the en- 
ergy density of the fluids are pi = E i /a 3 ( 1+ai > , where the 
constant of motion Ei is proportional to the total energy 
of the respective fluid. The corrections in Friedmann's 
equations resulting from the bulk-brane interaction are 
quadratic in the energy-momentum tensor of the matter 
fields on the brane. The resulting dynamics is noninte- 
grable and chaotic if the mass of the inflaton m ^ 0, 
allowing for metastable configurations that realize infla- 
tion due to parametric nonlinear of KAM tori present in 
the phase space of the model. 

Our analysis is restricted to dynamical configurations 
in which the dynamics of the scale factor is initially 
bounded in a potential well arising in the gravitational 
sector (a, p a ) due to the effective cosmological constant, 
the positive spatial curvature of the brane, and the bulk- 
brane corrections (connected with the perfect fluid com- 
ponents) that act as infinite potential barrier and is re- 
sponsible for the avoidance of the singularity a = 0. 
These configurations have the theoretical advantage over 
one single bounce models in that they avoid the prob- 
lem of initial conditions for the universe at past infin- 
ity. Furthermore they are past-eternal, oscillating about 
a stable Einstein universe configuration, that has no 
classical analogue and is favoured by maximum entropy 
considerations 15|. The skeleton of the phase space dy- 
namics is analyzed through its basic structures as critical 
points (stable and unstable Einstein universe), invariant 
plane, separatrices and deSitter attractors at infinity, as 
well as the foliation by KAM tori of the phase space in 
the nonlinear neighborhood of the stable Einstein uni- 
verse. The stable Einstein universe is a critical point of 
the dynamics when sourced by perfect fluids, or a limiting 
one-dimensional torus (namely, the topological product 
of a point times a minimal periodic orbit in the scalar 
field sector) when sourced by a scalar field in pure scalar 
field cosmologies. 

For particular domains of the parameter space 
(A4, <7, Ei , m) of the model, denoted windows of res- 
onance, these oscillatory bounded configurations turn 
metastable due to nonlinear parametric resonance phe- 
nomena, allowing the emergence of the universe into an 
inflationary phase. For numerical/analytical simplicity, 
we made an extensive examination of parametric reso- 
nance phenomena in braneworld models restricted to a 
dark energy component and radiation (i = dust,rad), 
plus the massive inflaton field and a dark energy com- 
ponent described by the effective cosmological constant 
in the brane. We constructed the resonant chart for the 
case of -Edust = -Erad = 10~ 3 , shown in Fig. 3, where 
the windows of resonance are the gray sheets. The white 
regions correspond to stable motion. When the system is 
driven towards a resonance window, nonlinear resonance 
of KAM tori takes place resulting in a complex dynamics. 
KAM tori that trapped the orbits are disrupted, and the 
initially bounded, oscillatory orbits of the system may 
escape to the deSitter attractor at infinity realizing in- 
flation. We illustrated this behaviour by Poincare maps 



in the plane (<£>',(/?), the origin of which corresponds to 
a stable or unstable periodic orbit whether respectively 
the system is in a stable region or a resonance window of 
the resonance chart. 

Each resonance window is characterized by an integer 
n > 2 and its main feature is the bifurcation of the sta- 
ble periodic orbit at the origin {ip — 0,p v — 0) into an 
unstable periodic orbit accompanied by one or two char- 
acteristic stable periodic orbits whether respectively n is 
odd or even. As the initial conditions of the expectation 
values ip are assumed to be small, being taken near the 
invariant plane ip = 0, p^ = 0, it follows that the para- 
metric domains of resonance are the ones that allow for 
inflation in the system. In this sense, since inflation is 
a sound paradigm for cosmology strongly sustained by 
observations, and if our present Universe is actually a 
braneworld, then the values of the cosmological parame- 
ters must be constrained to the resonance windows - with 
the braneworld inflated from initial conditions connected 
to a particular resonance. In particular, for fixed m and 
Ei, the brane tension a that regulates the strength of 
the effective Newton's constant in the brane will be re- 
stricted to small sheets depending on the integer n > 2 
(cf. Fig. 3). Therefore the larger the order of the reso- 
nance the stronger the gravitational coupling strength in 
the in the respective brane inflated due to a specific reso- 
nance. In this instance, we observe a quantization of the 
brane tension and consequently of the effective Newton's 
constant. 

The volume of the resonance windows are small as com- 
pared to the whole volume of the parameter space. The 
size of these 4-dim volumes depends strongly on the fluid 
content of the model. In Fig. 5 we exhibited the win- 
dows corresponding to the n = 3 resonance for pure dust, 
pure radiation and dust plus radiation in equal amounts 
(-E'rad = -Edust)- The case in which the perfect fluid com- 
ponent is dust only presents resonance windows with a 
reduced volume compared with pure radiation windows, 
analogously to the substructure of disruptive resonances 
in each window. This suggests that if we demand that 
our preinflationary models contain a component of dark 
matter in the form of dust this amount should be properly 
bounded in order that the model could realize inflation. 

Three distinct dynamical patterns are set up by the 
resonance, according to substructures in the resonance 
windows in models with -Edust < E ra d as observed by a 
detailed numerical examination. If the initial conditions 
correspond to configurations in the left border of the res- 
onance windows we have short time disruption of the 
initially bounded orbit with a rapid escape to inflation. 
On the other hand, if the initial conditions correspond 
to configurations in an intermediate threshold region of 
the window the orbit undergoes a long time diffusion in 
the stochastic sea surrounding the two secondary KAM 
stability islands, with posterior escape to inflation. Be- 
yond the threshold region, near the right border of the 
window, the orbit undergoes diffusion through large re- 
gions of phase space but remains bounded for times larger 
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than t = 100,000. This latter region of the resonance 
window should then be excluded as physically not ad- 
missable since they correspond to configurations that do 
not realize inflation. 

We also examined cosmological scenarios with a pure 
scalar field. Due to the absence of any perfect fluid com- 
ponent, the bulk-brane correction term that allows the 
dynamics to avoid the singularity is crucially dependent 
on the nonzero initial amplitude of inflaton through the 
conserved energy A = (^'(O) 2 + <p(0) 2 )/2. Contrary to 
the cases with a perfect fluid component, the stable Ein- 
stein universe is not a critical point but a limiting 1-dim 
tori, sourced by a scalar field. For a fixed a each en- 
ergy surface contains just one torus generated by a sin- 
gle orbit. The escape to inflation in the nonintegrable 
case is not associated to a parametric resonace pattern 
and occurs just for configurations of the scalar field with 
large initial conditions, in the small stochastic sea on the 
border of the main KAM island shown in Fig. 9. In 
this outer border we observe a mechanism of dynamical 
partial confinement of orbits, leading to finite time os- 
cillations before escape to inflation. Also, contrary to 
models with a perfect fluid component, the structure of 
the bouncing dynamics is extremely sensitive to the ini- 
tial amplitude and to the mass of the scalar field, and 
the presence of dynamical potential barriers allowing for 
bounces in the scale factor appears as a new feature of 
the dynamics. 

In Section VIII we examined the chaotic exit to in- 
flation for initial condition sets (corresponding to ini- 
tially expanding universes) taken in a small neighbor- 
hood about the stable separatix S. These sets are shown 
to have fractal basin boundaries connected to the code 
recollapse/escape to inflation leading to a chaotic exit to 



inflation. The fractality of the initial condition sets in- 
creases with to, namely, with the nonintegrability and 
their escape to inflation takes place smoothly involving 
no strong amplification of the inflaton field, as opposed to 
the metastable nonlinear resonance regime immediately 
before the exit to inflation, as discussed in Sections V-VI. 
We also observe the phenomenum of draining of these ini- 
tial condition basins from recollapse to escape behaviour, 
as time increases. For r — > oo only the homoclinic inter- 
section manifold remains in recurrent oscillatory motion. 
The homoclinic chaotic exit to inflation appears to be 
enhanced in the case of dust plus radiation, relative to 
the cases with pure dust or pure radiation components, 
contary to what occurs in the parametric resonance es- 
cape to inflation. We have not examined the possibility 
of chaotic escape to inflation in pure scalar field cosmolo- 
gies. 

Typically variation of the parameters can shrink or 
stretch the resonance zones. However the underlying 
pattern of resonance windows is maintained as we have 
checked numerically. In this sense the pattern is said to 
be structurally stable. Finally if our actual Universe is 
a brane inflated by a parametric resonance mechanism 
triggered by the inflaton, some observable cosmological 
parameters (e.g. Newton's gravitational constant on the 
brane and the mass of the inflaton) should then have 
a signature of the particular resonance from which the 
brane inflated and consequently of the particular value 
of the parameters [a, to, -Edust, -Brad) that were favored in 
the early dynamical regime of the universe. 
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